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We report on the calculation of the next-to-leading order QCD corrections to the production of ti 
pairs in association with two hard jets at the Fermilab TeVatron and CERN Large Hadron Collider. 
Results for the integrated and differential cross sections are given. The corrections with respect to 
leading order are negative and moderate. A study of the scale dependence of our NLO predictions 
indicates that the residual theoretical uncertainty, due to higher order corrections, is 21% for the 
TeVatron and 15% for the LHC. In case of the TeVatron, the forward-backward asymmetry of the 
top quark is calculated for the first time at next-to-leading order. With the inclusive selection of 
cuts, this asymmetry amounts to -4 fb ,lo = —10.3% at leading order and -4fb,nlo = —4.6% at 
next-to-leading order. All results presented in this paper have been obtained with the help of the 
Helac-Nlo package. 

PACS numbers: 12.38.Bx, 14.65. Ha, 14.80.Bn 



in: 
43 



in 

00 
(N 

od 
o 



X 



I. INTRODUCTION 

Top quark production in association with two jets con- 
stitutes an important background for many new parti- 
cle searches at the Tevatron and at the LHC. Promi- 
nent examples include searches for Higgs boson decays 
H — > WW* and H — > bb, where the Higgs boson is pro- 
duced via weak boson fusion or via an associated produc- 
tion with a ti pair respectively. Higgs boson production 
in association with top quarks has a strongly decreasing 
cross section with increasing Higgs boson mass, which 
makes the process useful only in the low mass range, 
mg < 135 GeV, when Higgs boson decays in bb are im- 
portant. This potential discovery channel gives a unique 
access to the top and bottom Yukawa couplings. Com- 
bined analyzes of ATLAS [J and CMS gf have shown 
that a 3cr evidence of the signal above the dominant ttbb 
and tijj backgrounds can be obtained for Mh < 130 
GeV at the LHC, if enough luminosity, C = 60 fb _1 , is 
collected. However, a reconstruction of the H — > bb mass 
peak is difficult because of an identification problem in 
the signal. The bb pair can be chosen incorrectly, due to 
a lack of a distinctive kinematic feature of jets from the 
Higgs boson decay. Moreover, the b-tagging efficiency 
takes a crucial part in this kind of analyzes, since two b- 
jets for a Higgs boson candidate can arise from mistagged 
light jets. Therefore, a very precise knowledge of the 
backgrounds is necessary, if the tiH channel is to be of 
any use. 

On the other hand, a detailed ATLAS examination 
of the vector boson fusion channel at the LHC, pp — > 
Hjj — > W + W~jj — > Ivtvjj, has shown that a signifi- 
cance larger than 3(7 can be obtained for a luminosity of 
C = 10 fb _1 in the Higgs boson mass range above ttih ~ 
130 GeV, i.e. when BK(H -> WW*) is large enough 0. 
When combining this channel with the Ivjj mode, one 



can even obtain a ~ 5ct significance for the above mass 
range. In fact, the pp — > Hjj — > WW jj channel becomes 
very powerful at higher Higgs boson masses. One can ar- 
rive at a signal significance that is above 5er in a wide 
mass range, namely from 140 GeV up to 190 GeV, for 
the same luminosity. Although feasible and competitive, 
this channel might prove to be rather difficult since the 
Higgs boson mass peak cannot be directly reconstructed. 
Consequently, one cannot measure background processes 
from the side bands. The most important background 
processes, ttjj and QCD W + W~jj, need to be known, 
therefore, as precisely as possible. 

Meanwhile at the TeVatron, with up to 7.1 fb _1 of data 
analyzed at CDF, and up to 8.2 fb _1 at DO, the Higgs 
boson mass range of 158 GeV < rnn < 173 GeV has 
been excluded [4|. Very recently, the LHC experiments 
have presented their Higgs boson search results based on 
1 fb _1 of analyzed data. CMS has excluded the Stan- 
dard Model Higgs boson in the 149 GeV - 206 GeV and 
300 GeV - 440 GeV windows 0, while ATLAS has ex- 
cluded the 155 GeV - 190 GeV and 295 GeV - 450 GeV 
windows [f|. The low mass exclusion is dominated by 
the search of the H — > W + W~ — > Ivlv final state, while 
the high mass one is dominated by H — > ZZ with three 
different combinations of Z decays. Considering these re- 
cent experimental results, an excellent understanding of 
lighter Higgs boson search channels is more timely than 
ever. 

Apart from its significance as a background process 
to various new physics searches it turns out that ttjj 
production can also be an important signal process. With 
the total integrated luminosity accumulated so far at the 
TeVatron around 3000 events should have been collected 
so far by each experiment. As for the LHC, in June of this 
year both experiments ATLAS and CMS reached 1 fb _1 
of collected data. This can be directly translated to about 
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10000 events already stored. Since they are in the process 
of collecting data, the tijj process can soon become of 
great importance for a more precise understanding of the 
top quark events topology. 

However, it is a well known fact that for processes in- 
volving strongly interacting particles, the leading order 
(LO) cross sections are affected by large uncertainties 
arising from higher order corrections. If, at least, next- 
to-leading order (NLO) QCD corrections to these pro- 
cesses are included, total cross sections can be denned in 
a reliable way. The NLO QCD corrections to pp — > ttbb 
production have already been calculated by two inde- 
pendent groups 043; and perfect agreement has been 
found. In addition, NLO corrections to the remai ning 
two background processes, namely pp — > W + W~jj [10| 
and pp — > tijj |ll| have recently been evaluated. In 
this paper, we have extend our previous study on the 
NLO QCD corrections to the pp — > tijj production at the 
LHC, where the integrated cross section together with the 
scale dependence of the total cross section for one partic- 
ular set of cuts for the center of mass system energy of 
y/s = 14 TeV have been studied. Moreover, we have pre- 
sented there three differential distributions, the invariant 
mass distribution of the first and the second hardest jet 
together with the transverse momentum distributions of 
the first and the second hard jet. We supplement this 
discussion here with more results. In particular, results 
for the TeVatron case are going to be presented for the 
first time. In case of the LHC, we use the current center 
of mass system energy of y/s = 7 GeV. Additionally to 
the integrated cross sections for different setups, various 
differential distributions are given. 

The article is organized as follows. In Section [TT] we 
briefly summarize the methods we use to obtain our NLO 
QCD results. In Section [Hi] numerical results are pre- 
sented, first for the TeVatron and later on for the LHC 
case. We give our conclusions in Section IIVI Finally, in 
an Appendix we give values of squared matrix elements 
for the virtual corrections for one phase space point. 



II. FRAMEWORK OF THE CALCULATION 

The calculation of NLO corrections to hadronic top 
pair production in association with two jets proceeds 
along the same lines as our earlier work on pp — > tibb 
II, pp -> ttjjjX^ and pp(pp) tt -> W+W-bb -> 
e + v e fi~v f _ l bb [lJI- The methods developed there have 
therefore been straightforwardly adapted and only need 
a brief recollection here. Let us emphasize, that all parts 
of the calculation are performed fully numerically in a 
completely automatic manner. 

A. Tree level contributions 

At tree level the it production in association with 
two jets proceeds via the scattering of two gluons, two 



TABLE I. Partonic subprocesses contributing to the lead- 
ing order process pp(pp) — > tijj at 0{oQ). The number of 
Feynman diagrams corresponding to these subprocesses is also 
shown. 



Partonic 
Subprocess 


Number Of Feynman 
Diagrams 


99 -> tigg 


123 


99 ttqq 


36 


qq -i tigg 


36 


91 -» ttqg 


36 


qg -> ttqg 


36 


qq' — > tiqq' 


7 


qq — > tiqq 


7 


qq — > tiqq 


14 




FIG. 1. A representative set of Feynman diagrams contribut- 
ing to the LO hadronic tijj production at O(a^). Double 
brown lines correspond to top quarks, single lines to light 
quarks and wiggly ones to gluons. 



quarks or quark and gluon. All contributions of the or- 
der of 0{a^) can be grouped into eight easily identifiable 
topologies presented in the Table [TJ Both q and q 1 span all 
quarks and anti-quarks. Also shown are the numbers of 
Feynman diagrams corresponding to these subprocesses. 
We emphasis here that we do not use Feynman diagrams 
in our calculations. Instead, matrix elements and total 
cross sections are obtained with an iterative algorithm 
based on the Dyson- Schwinger equations [l3|, [lij]. Never- 
theless, we present the number of Feynman diagrams as a 
measure of complexity. All partonic channels have been 
calculated with the Helac-Dipoles HH package. An 
independent calculation using the Helac-Phegas [iH - 
[l8j event generator has been performed. Let us stress 
that Helac-Phegas has been already extensively used 
and tested, see e.g. [l9l - l24T |. Moreover, both libraries 
compute relevant contributions in a very efficient way, 
using off-shell recursive equations and avoiding a multi- 
ple evaluation of recurring building blocks. In every case 
we have found full numerical agreement. The optimiza- 
tion and phase space integration was executed with the 
help of Parni [2f| and Kaleu [2(J. A few examples of 
LO graphs are shown in Figure [TJ 



TABLE II. The number of one-loop Feynman diagrams for 
the pp(pp) — > tijj process at 0{a 5 s ). 





IN U lvli3il/rt yjr r Fj i IN 1V1A1N 

T^T APT! A 1V/TQ 
LJ 1 A^j xt J\ IV! O 


- 

gg —¥ ttgg 


401U 


gg -t tiqq 


1100 


qq -t ttgg 


1100 


gq -> tiqg 


1100 


qg -> tiqg 


1100 


qq' —> tiqq' 


205 


qq -> tiq'q' 


205 


qq — > tiqq 


410 



B. Virtual corrections 

The virtual corrections are obtained from the interfer- 
ence of the sum of all one-loop diagrams with the Born 
amplitude. One can classify them into self-energy, vertex, 
box-type, pentagon-type and hexagon-type corrections. 
In Table HI the number of one-loop Feynman diagrams 
for the pp(pp) —> tijj process at 0{a b s ) as obtained with 
QGRAPH [27| is given. Again they are shown as a measure 
of complexity of our calculation since we are free from the 
task of computing Feynman graphs for a given process. 
Typical examples of virtual graphs are shown in Figure 
[21 The virtual corrections are calculated in d = 4 — 2e 
dimensions in the 't Hooft-Veltman [28[ version of the di- 
mensional regularization scheme. The singularities stem- 
ming from infrared divergent pieces are canceled by the 
respective poles in the integrated counter terms of the 
dipole subtraction approach. The finite contributions of 
the loop diagrams are evaluated numerically in d = 4 di- 
mension. We monitor the numerical stability by checking 
Ward identities at every phase space point. The events 
which violate gauge invariance are not discarded from 
the calculation of the finite part, but rather recalculated 
with higher precision. In order to compute the relevant 
contributions the Helac-ILoop [29[ approach is used. 
It is based on the Helac-Phegas program to calculate 
all tree-order like ingredients and the OPP [3(| reduc- 
tion method. The cut-constructible part of the virtual 
amplitudes and of the rational terms, R\, are computed 
using the CutTools [3l[ code. The remaining part of 
the rational term, R2, is obtained by the use of extra 
Feynman rules as described in fliTi ■''>"> . The evaluation of 
scalar integrals is performed with the help of OneLOop 
[1^, [3(|. Renormalization is done as usual, by evaluating 
tree level diagrams with counterterms. For our process, 
we chose to renormalize the coupling in the MS scheme 
with five active flavors and the top quark decoupled, just 
as in [H, [ll[ , and the mass in the on-shell scheme (wave 
function renormalization is done in the on-shell scheme 
as it must be). 
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FIG. 2. A representative set of Feynman diagrams contribut- 
ing to the virtual corrections to hadronic tijj production at 
0(a 5 a ). Double brown lines correspond to top quarks, single 
lines to light quarks and wiggly ones to gluons. 



TABLE III. Partonic subprocesses contributing to the sub- 
tracted real emissions at 0(a s ) for the pp(pp) — > tijj process. 
Also, the number of Feynman diagrams and the number of 
dipoles corresponding to these subprocesses are presented. 



Partonic 


Number Of Feynman 


Number Of 


Subprocess 


Diagrams 


Dipoles 


99 -> ttggg 


1240 


75 


gg -> tiqqg 


341 


55 


qq -> tiggg 


341 


75 


gq -> tiq'q' q 


64 


25 


gq -> tiqgg 


341 


65 


qg -> tiq'q' q 


64 


25 


qg -> tiqgg 


341 


65 


qq' -> ttqq'g 


64 


40 


qq -t tiq'q' g 


64 


35 


qq -> tiqqg 


128 


45 



C. Real corrections 

The real emission corrections to the LO process arise 
from tree level amplitudes with one additional parton, 
an additional gluon, or a quark anti-quark pair replac- 
ing a gluon. All possible contributions of the order of 
O(o.g) can be divided into ten subprocesses which are 
given in Table IIII1 where again q and q' stand for up- 
or down- type quarks. Moreover, the number of Feyn- 
man diagrams and the number of dipoles corresponding 
to these partonic reactions are shown. We employ the 
dipole subtraction formalism in the form as proposed by 
Catani and Seymour [37| , to extract the soft and collinear 
infrared singularities and to combine them with the vir- 
tual corrections. Specifically, the formulation presented 
in [38| for massive quarks has been used with the ex- 
tension to arbitrary helicity eigenstates of the external 
partons [lj|, as implemented in Helac-Dipoles. After 
combining virtual and real corrections, singularities con- 
nected to collinear configurations in the final state as well 
as soft divergences in the initial and final states cancel 
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for collinear-safe observables automatically after apply- 
ing a jet algorithm. Singularities connected to collinear 
initial-state splittings are removed via factorization by 
parton distribution function redefinitions. In our case 
the cancellation of divergences has been checked numer- 
ically for a few phase space points. Let us mention here 
that similarly to the LO case, the phase space integra- 
tion is executed with the help of Kaleu |26[. However, 
in case of subtracted real emission a version with dipoles 
channels have been used instead, see [lI for details. 



III. NUMERICAL RESULTS 

In the following we present predictions for the pro- 
cess pp(pp) — > ttjj + X both at the TeVatron run II 
with yfs = 1.96 TeV and at the Large Hadron Collider 
(LHC) with y/s = 7 TeV. We have consistently employed 
the MSTW2008 [HI set of parton distribution functions 
(PDFs). In particular, we take MSTW2008LO PDFs 
with 1-loop running a s in LO and MSTW2008NLO PDFs 
with 2-loop running a s in NLO, including five active fla- 
vors. Contributions induced by bottom-quark density are 
taken into account. We use the Tevatron average mass 
of the top quark ra t — 173.3 GeV as measured by 
the CDF and DO experiments. The masses of all other 
quarks, including b quarks, are neglected. If not stated 
otherwise, the renormalization and factorization scales 
are set to ji r = jip = m t . Jets are defined via an in- 
frared safe algorithm, the kr (U-43], anti—kx Q (our 
default) and the inclusive Cambridge/ Aachen (C/A) al- 
gorithm [45| . In all cases no clustering can take place 
unless two partons of pseudorapidity 



TABLE IV. Integrated cross section at LO and NLO for 
pp — > ttjj + X production at the TeVatron run II with 
yfs = 1.96 TeV. Results for three different jet algorithms 
are presented, the anti-kr, kr and the Cambridge/ Aachen jet 
algorithms with R — 0.4. The scale choice is hr = /j,f = mt. 



Cuts 


o-lo [pb] 


^NLO [Pb] 


°nIo [Pb] ^nlo [Pb] 


p Tj > 20 GeV 








ARjj > 0.4 


0.3584(1) 


0.2709(5) 


0.2734(4) 0.2734(4) 


\Vi\ < 2.0 









TABLE V. Integrated cross section at LO and NLO for 
pp — > ttjj + X production at the TeVatron run II with 
y^s = 1.96 TeV. Results for three different jet algorithms 
are presented, the anti-kr , &t and the Cambridge/ Aachen jet 
algorithms with R — 0.8. The scale choice is fiR = (if = mt- 



Cuts 


o-lo [pb] 


anti — k-T r l i 
°NLO [Pb] 


^nlo [Pb] <%lo [Pb] 


p T] > 20 GeV 








ARjj > 0.8 


0.2876(1) 


0.2467(3) 


0.2494(3) 0.2491(3) 


\Vi\ <2.0 









AR j0 = 



A^. + Ayl > 0.4 



33 



(4) 



We will refer to this set of cuts as TeVatron default se- 
lection. All results presented in this Section are given 
in picobarns, except for the Table IVlTl where results are 
shown in femtobarns. 



\<1\ 



■In [tan (0/2)] < 5, 



(1) 



where 9 is the angle between the parton momentum and 
the beam axis, are separated by less than R in the (77, (f>) 
plane. Reconstructed jets are ordered in px- Their mo- 
menta are formed as the four- vector sum of massless par- 
tons. The outgoing (anti-)top quarks are left on-shell 
with unrestricted kinematics. They are assumed to be 
always tagged. 



1. Integrated cross sections 

Now we turn to a discussion of the numerical results. 
We start with the total cross sections for the TeVatron 
default selection. We find 



ill7(pp~>m = 0-3584 +°-«> pb , 



(5) 



TeVatron 



ttjj) 



27Q9 +0.0014(0.5%) , 
U - Z,Uy -0.0566(21%) t" J - 



(6) 



We calculate the partonic cross sections for events with 
at least two hard jets. Jets are obtained via the anti—kT 
jet algorithm with a resolution parameter R — 0.4 and 
are required to have 

PTj — \J Pxj + Pyj > 20 GeV, (2) 



TABLE VI. Integrated cross section at NLO for pp — > ttjj + 
X production at the TeVatron run II with yfs = 1.96 TeV. 
Results for two different values of the a ma x parameter, jet 
resolution R and ARjj cut for the anti-kr jet algorithm are 
presented. The scale choice is fj,R — [if = mt- 



Cuts 




a max =0.01 r t 1 
NLO [P^ 


a max = 1.00 r w 
NLO lP b J 




> 0.4, R = 0.4 


0.2709(5) 


0.2710(2) 




> 0.8, R = 0.8 


0.2467(3) 


0.2466(2) 
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TABLE VII. Integrated cross section at LO and NLO for 
pp — > tijj + X production at the TeVatron run II with 
\fs — 1.96 TeV. Results for the anti-kr jet algorithm with 
R — 0.4 are presented. In the last two columns the K, factor, 
defined as the ratio of the NLO cross section to the respective 
LO result, and NLO corrections in % are given. The scale 
choice is [ir = fiF = m-t- 



Cut 



olo [fb] 



anti — k-j 
a NLO 



[fb] 



Pt 3 > 20 GeV 
p Tj > 40 GeV 
p Tj > 60 GeV 
p Tj > 80 GeV 



358.4(1) 
79.29(3) 
23.85(1) 
8.274(3) 



271.2(4) 
54.11(9) 
14.55(3) 
4.520(13) 



0.76 
0.68 
0.61 
0.55 



-24 
-32 
-39 
-45 





0.25 r 




0.2 


On 


0.15 






c 


0.1 


X3 





0.05 




-1.5 -1 



At the central value of the scale, fj,R — /if = Tnti the 
full cross section receives moderate NLO correction of 
the order of —24%. The scale dependence is indicated 
by the upper and lower indices. The upper (lower) index 
represents the change when the scale is shifted towards 
fx = mt/2 {(i = 2m t ). Rescaling the common scale from 
the default value rxit up and down by a factor 2 changes 
the cross section at LO by about 94%. Through the 
inclusion of NLO QCD corrections the scale uncertainty 
is dramatically reduced down to 21%. 

To assess the effect of changing the jet algorithm, we 
compare NLO results for the hr, anti—kx and the in- 
clusive Cambridge/ Aachen (C/A) jet algorithms as pre- 
sented in Table ITVl With the jet resolution parameter of 
R = 0.4 the difference is within 1%. Total cross sections 
with a higher jets separation cut, ARjj > 0.8, and the jet 
resolution parameter R — 0.8 are shown in Table IVl The 
NLO QCD corrections are reduced down to —16% in that 
case. This is mostly due to a large, —20% change, in the 
LO cross section as compared to a rather moderate, —9% 
shift at NLO. With an increase of the R parameter the 
perturbative difference between the algorithms should be 
more visible because the phase space for inclusion of an 
extra parton in the jet is much larger. Nevertheless in 
our case it remains within 1%. 

For completeness and as part of checks, the integrated 
NLO cross sections for two different values of the a max 
parameter are given in Table I VII Let us remind the 
reader that the a max parameter is a common modifica- 
tion of subtraction terms in the phase space away from 
the singularity. We consider two extreme choices, namely 
q m p = 1, which corresponds to the original formulation 
of [38|, and a max = 0.01. As can be observed, the in- 
dependence of the final result on the value of a max is 
obtained at the permil level. 

In Table IVIII we provide LO and NLO predictions 
for the integrated cross sections for the TeVatron de- 
fault selection, however, for different values of the trans- 
verse momentum cut for the first and the second hard- 
est jet. The values are given for the central scale value 
I^R = (J'F = TUf Up to ~ 50 GeV, rather moderate cor- 




-1.5 -1 -0.5 0.5 1 1.5 

y t 

FIG. 3. Differential cross section distributions as a function 
of rapidity, y t , of the top and anti-top quark at LO (upper 
panel) and NLO (lower panel) for pp — > tijj + X production 
at the TeVatron run II with *Js = 1.96 TeV. The (orange) 
solid curve corresponds to the top quark, whereas the (brown) 
dashed one to the anti-top quark. 



rections are acquired, however, for the high value of pt 
cut, a substantial increase is noticed. For apTj change by 
a factor of 4, NLO QCD corrections are almost doubled. 
Compared to the total NLO it cross section, 

a^oiPP^ti) = 6.68+°; 3 7 ^ o) pb, (7) 

obtained for the same setup, we find that for the small 
PTj cut of 20 GeV the tijj events represent only 4% of the 
total cross section. The fraction is reduced substantially, 
down to 0.8%, 0.2% and 0.07% respectively, if a px 3 cut 
of 40 GeV, 60 GeV and 80 GeV is chosen instead. 



2. Forward- Backward Asymmetry 

In the next step we calculate the top quark forward- 
backward asymmetry. At LO the asymmetry is defined 
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-1.5 -1 -0.5 0.5 1 1.5 

y t 




-20 



-40 

-1.5 -1 -0.5 0.5 1 1.5 

y t 

FIG. 4. Differential charge asymmetry, A(yt), as a function 
of the (anti-)top quark rapidity at LO (upper panel) and NLO 
(lower panel) for pp — > tijj + X production at the TeVatron 
run II with = 1.96 TeV. NLO\ refers to a result with 
a consistent expansion in a s , while NLO2 to the unexpanded 
one. 




At _ a hO ( 1 i ^°"nLO ^°"nLO ^ f-\ n \ 

-^fb,nlo — ~ r~ 1 ^ = + — • \ w ) 

ct lo V <y w <j lo J 

The size of the forward-backward asymmetry of the top 
quark at LO for tijj production amounts to 

4' — — fl 1 +0 003 /1 1\ 

•^FB^O — u - lud -0.004- V iJ J 

The NLO correction to the asymmetry as calculated from 
Eq. (|9]), and with a consistent expansion in a s reads 

■4b, N lo = -0.046 ±S:8Si - (12) 

Had we used an unexpanded ratio of the NLO cross sec- 
tions, the result would rather be 

-4fb,nlo = ~0.0581q'o42 j (13) 

which means that the two definitions give reasonably con- 
sistent results for the central scale. However, the theo- 
retical error as calculated from the scale dependence is 
almost as large as the asymmetry itself in the latter case. 

In Figure [3] rapidity distributions for the top and anti- 
top quarks are presented at LO and NLO. In both cases 
results are not symmetric around yt = and are shifted 
to a larger yt for the anti-top quarks and smaller yt for the 
top quarks. This shows that anti-top quarks are prefer- 
entially emitted in the direction of the incoming protons 
(forward direction by definition), as implied by the sign 
in Eq. [T2] In Figure [H we have also plotted the differ- 
ential charge asymmetry, A(yt)- In the NLO case both 
results with and without a consistent expansion in a s 
are shown. In both cases, it is clearly visible that the 
forward-backward asymmetry of top quarks in tijj pro- 
duction is significantly reduced if NLO QCD corrections 
are taken into account. 



as 

A t _ Vhojyt > 0) - cr LO (y t < 0) 

-^FB.LO — / ~ n \ I / T7TT ' \°J 

0x0(2/* > 0) + a w (yt < 0) 

where = ctlo(j/* > 0)±<7Lo(z/t < 0) is evaluated with 
LO PDFs and LO a s . On the other hand, the asymmetry 
at NLO is defined by 

At _ a LO + ^°"nLO (q\ 

-^fb.nlo — ~f ; ; + > W 

CT LO + ^NLO 

where Sa^ LO are the NLO contributions to the cross sec- 
tions and crj~ are eva l ua ted with NLO PDFs and NLO 
a s . In case of large differences between the LO and NLO 
asymmetries, it is necessary to expand Eq. ^ to first or- 
der in a s . Indeed, the ratio in Eq. (|9|) generates contribu- 
tions of 0{a 2 s ) and higher, which are affected by the un- 
known next-to-next-to-leading order contributions. The 
following definition has been used, for the same reasons, 



3. Differential cross sections 

In addition to the normalization of the integrated cross 
section, radiative corrections can affect the shape of vari- 
ous kinematic distributions. To quantify the size of these 
distortions we have tested a variety of differential distri- 
butions. 

We start in Figure [5] with the differential distribution 
in the rapidity of the top and anti-top, yt, and in Figure 
[5] with an averaged transverse momentum of the top and 
anti-top quarks. The dash-dotted curve corresponds to 
the LO, whereas the solid one to the NLO result. The 
upper panels show the distributions themselves and ad- 
ditionally include the scale-dependence bands obtained 
with a scale variation by a factor of two. The lower pan- 
els display the ratio of the NLO value to the LO result, 
for the central scale values of /i = (1r = /ip = m t , called 
the differential K, factor. For the angular distributions we 
observe negative corrections of the order of 15% — 30%. 
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-1.5 -1 -0.5 0.5 1 1.5 -1.5 -1 -0.5 0.5 1 1.5 

y t y- t 

FIG. 5. Differential cross section distributions as a function of the rapidity of the top (left panel) and anti-top (right panel) for 
pp — > tijj + X production at the TeVatron run II with y/s = 1.96 TeV. The dash-dotted curve corresponds to the LO, whereas 
the solid one to the NLO result. The uncertainty bands depict scale variation. The lower panels display the differential fC factor. 
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FIG. 6. Differential cross section distributions as a func- 
tion of the averaged transverse momentum of the (anti-)top 
for pp — >■ tijj + X production at the TeVatron run II with 
\fs — 1.96 TeV. The dash-dotted curve corresponds to the 
LO, whereas the solid one to the NLO result. The uncer- 
tainty bands depict scale variation. The lower panels display 
the differential K, factor. 



More precisely, 15% corrections are observed in the for- 
ward (backward) regions of a top (anti-top) distribution 
while 30% corrections are in the backward (forward) re- 
gions. For the transverse momentum distribution up to 



100 GeV moderate and negative corrections of the order 
of 20% are reached while the tail exhibits even —60% 
corrections. 

In Figure [7] and Figure [5] differential cross section dis- 
tributions as a function of rapidity, y t t, transverse mo- 
mentum, PT tt , and invariant mass, m t i, of the ti pair are 
presented. Also shown is the differential distribution in 
the total transverse energy of the system, Ht- The latter 
is defined as 

H T = m Tt + m T{ + ^ pr h , (14) 

i=l,2 

where 



m Tt = \Jm 2 t +p Tt . (15) 

The rapidity distribution of the ti pair has almost con- 
stant, negative and moderate 20% — 25% corrections. 
Both, m t t andpT tt - distributions get large negative 50% — 
60% corrections for the higher values of these observ- 
ables, which have to be compared with also negative, 
but 20% — 25% corrections close to the ti threshold for 
m t i and for smaller values of the transverse momentum 
of the ti pair for the px t{ distribution. At the begin of the 
Ht spectrum negative and small 10% — 15% corrections 
can be seen. However, at the same time even —70% can 
be reached for higher values. 

In the next step, the jet kinematics is presented. We 
start with transverse momentum and rapidity distribu- 
tions of the 1 st and the 2 nd hardest jet as depicted in 
Figure [D and Figure [TUJ Again, in case of rapidity dis- 
tributions, negative and almost constant 20% — 25% cor- 
rections are obtained. On the other hand, also negative 



8 




% 

o 



0.008 



| 0.006 




50 100 150 200 250 
P T -[GeV] 



FIG. 7. Differential cross section distributions as a function of rapidity (left panel) and transverse momentum (right panel) of 
the tt pair for pp — » tijj + X production at the TeVatron run II with yfs — 1.96 TeV. The dash-dotted curve corresponds to 
the LO, whereas the solid one to the NLO result. The uncertainty bands depict scale variation. The lower panels display the 
differential K, factor. 
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FIG. 8. Differential cross section distributions as a function of the invariant mass of the tt pair (left panel) and of the total 
transverse energy (right panel) for pp — > ttjj + X production at the TeVatron run II with t/s — 1.96 TeV. The dash-dotted 
curve corresponds to the LO, whereas the solid one to the NLO result. The uncertainty bands depict scale variation. The lower 
panels display the differential K, factor. 



40%— 50% corrections are visibled at the tails oipTj spec- 
trums, as compared to minus 20% in the low transverse 
momentum range. 

And finally, in Figure [TTl the invariant mass, rrijj, of 
the 1 st and the 2 nd hardest jet together with their sep- 
aration in the (y, <p) plane are shown. In both cases al- 
most constant and rather moderate differential /C-factors 



of 0.75 can be seen. Except for the tails of Ai?^- dis- 
tribution where distortions are much higher. However, 
fluctuations visible there arc a reflection of the limited 
statistics of the Monte Carlo integration. At the TeVa- 
tron, ARjj is dominated by the transverse plane peak 
at the begin of the spectrum, which suggests that two 
hardest jets are preferably produced almost collinear to 
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FIG. 9. Differential cross section distributions as a function of transverse momentum (left panel) and rapidity (right panel) of 
the 1st hardest jet forpp — > tijj + X production at the TeVatron run II with y/s = 1.96 TeV. The dash-dotted curve corresponds 
to the LO, whereas the solid one to the NLO result. The uncertainty bands depict scale variation. The lower panels display the 
differential K, factor. 
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FIG. 10. Differential cross section distributions as a function of transverse momentum (left panel) and rapidity (right panel) of 
the 2nd hardest jet for pp — > tijj + X production at the TeVatron run II with yfs = 1.96 TeV. The dash-dotted curve corresponds 
to the LO, whereas the solid one to the NLO result. The uncertainty bands depict scale variation. The lower panels display the 
differential K, factor. 



each other. 

Overall, we can say that at the TeVatron employing a 
fixed scale /i = = fiF = mti the NLO corrections to 
transverse momentum distributions are rather substan- 
tial. Moreover, they do not simply rescale the LO shapes, 
but induce distortions at the level of 40%. The same ap- 
plies to the invariant mass distribution of the tt pair. In 



case of the total transverse energy, even 50% — 60% de- 
formations can be observed. For angular distributions 
and the invariant mass of two tagging jets we observe 
negative and rather moderate corrections of the order of 
15% — 30%, which turn out to be relatively constant. 
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FIG. 11. Differential cross section distributions as a function of the invariant mass of the 1st and the 2nd hardest jet (left panel) 
and ARjj separation (right panel) for pp — > tijj + X production at the TeVatron run II with yfs = 1.96 TeV. The dash-dotted 
curve corresponds to the LO, whereas the solid one to the NLO result. The uncertainty hands depict scale variation. The lower 
panels display the differential K, factor. 



TABLE VIII. Integrated cross sections at LO and NLO for 
pp — > tijj + X production at the LHC with y/s = 7 TeV. 
Results for three different jet algorithms are presented, the 
anti-kr, kr and the Cambridge/Aachen jet algorithms with 
R = 0.5. The scale choice is /j,r — fiF = mt. 

Cuts <tlo [pb] o" NLO 

kT [Pb] 

°NLO 

[Pb] °"nlo [Pb] 



p Tj > 50 GeV 

ARjj > 0.5 13.398(4) 9.82(2) 9.86(2) 9.86(2) 
ll/j|<2.5 



TABLE IX. Integrated cross section at LO and NLO forpp — > 
tijj + X production at the LHC with y/s — 7 TeV. Results 
for three different jet algorithms are presented, the anti-kr , 
kr and the Cambridge/Aachen jet algorithms with R = 1.0. 
The scale choice is [ir — fiF — mt- 

Cuts cr LO [pb] o-j* r n L o" kT [pb] 0& o [pb] a°[% [pb] 



p Tj > 50 GeV 

AR j:i > 1.0 11.561(4) 9.95(2) 10.06(2) 10.04(2) 
M <2.5 



B. LHC 

As in the TeVatron case we calculate the partonic cross 
sections for events with at least two hard jets. Jets are 
obtained via the anti—kr jet algorithm with a resolution 
parameter R — 0.5 and arc required to have 

PTj > 50 GeV , \ yj \ < 2.5 , ARjj > 0.5 . (16) 

We will refer to this set of cuts as LHC default selec- 
tion. All results presented in this Section are given in 
picobarns. 

1. Integrated cross sections 

We turn now our attention to the LHC results. Let us 
stress here, that the cross sections and differential distri- 
butions presented in this section do not coincide with the 



TABLE X. Integrated cross section at NLO for pp — > tijj + 
X production at the LHC with yfs = 7 TeV. Results for 
two different values of the a m ax parameter and ARjj cut are 
presented. They are obtained with the anti-kr jet algorithm. 
The scale choice is (ir = (j,f = nit- 



Cuts 




Q, m , x =0.01 r U l 
^NLO [Pbj 


°NLO [Pbj 


AR J3 


>0.5, R = 0.5 


9.82(2) 


9.81(1) 


ARjj 


> 1.0, R = 1.0 


9.95(2) 


9.94(1) 



ones shown in [11]. In the latter publication, a different 
center of mass energy, yfs — 14 TeV, has been assumed. 

First, integrated cross sections with the LHC default 
selection as defined in the previous section, are provided. 
We obtain 

o-t* c (pp -+ tijj) = 13.398 tl^lllS pb, (17) 
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FIG. 12. Differential cross section distributions as a function of the rapidity of the top quark (left panel) and of the anti-top 
quark (right panel) for pp — > tijj + X production at the LHC with •Js = 7 TeV. The dash-dotted curve corresponds to the LO, 
whereas the solid one to the NLO result. The uncertainty bands depict scale variation. The lower panels display the differential 
K, factor. 



TABLE XI. Integrated cross section at LO and NLO for pp — > 
tijj -\-X production at the LHC with yfs = 7 TeV. Results for 
the anti-kr jet algorithm with R — 0.5 are presented. In the 
last two columns the K, factor, defined as the ratio of the NLO 
cross section to the respective LO result, and NLO corrections 
in % are given. The scale choice is fin — [if = mt- 



where the central value refers to fiR — [If = mt, the up- 
per value to fi — m t /2 and the lower value to /i = 2mt- 
At the central value of the scale the total cross section 
receives negative and moderate corrections of the order 
of 27%. Similar to the TeVatron case, the improvement 
in the scale stability at the NLO is prominent. The NLO 
cross section is very stable against changes in the renor- 
malization and factorization scales, in contrary to the LO 
result where an 87% scale dependence is obtained. 

As in the TeVatron case we have checked the effect 
coming from the jet algorithm variation. Results are pre- 
sented in Table IVIIII For the jet resolution parameter 
R = 0.5 the difference between hp, anti — kj< and the 
inclusive Cambridge/ Aachen (C/A) jet finders is within 
0.4%. With an increase of R values to R = 1.0 the differ- 
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$ 0.075 

Q. 





Cut 


o"lo [pb] 


%lo [P b ] 


K 


[%] 


H 
Q. 

b 


0.05 


PTj 


> 50 GeV 


13.398(4) 


9.82(2) 


0.73 


-27 




0.025 




> 75 GeV 


5.944(2) 


4.115(8) 


0.69 


-31 






PTj 


> 100 GeV 


3.018(1) 


1.944(4) 


0.64 


-36 







PTj 


> 125 GeV 


1.665(1) 


0.993(2) 


0.60 


-40 




1.5 




LHC 
CT NLO 


{PP -> tijj) = 


q S9 -1.48(15%) 
y ' 0Z -1.47(15%) 


pb. 


(18) 
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FIG. 13. Differential cross section distributions as a function 
of the averaged transverse momentum of the top and anti- 
top for pp — > tijj + X production at the LHC with ^fs = 
7 TeV. The dash-dotted curve corresponds to the LO, whereas 
the solid one to the NLO result. The uncertainty bands depict 
scale variation. The lower panels display the differential K, 
factor. 



ence moves upwards to 1.1% as can be seen from Table 



IIXI Moreover, for a higher jet separation cut ARjj > 1.0 
and a higher jet resolution parameter R = 1, NLO QCD 
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FIG. 14. Differential cross section distributions as a function of rapidity (left panel) and transverse momentum (right panel) of 
the ti pair for pp — > tijj + X production at the LHC with ^/s — 7 TeV. The dash-dotted curve corresponds to the LO, whereas 
the solid one to the NLO result. The uncertainty bands depict scale variation. The lower panels display the differential K, factor. 
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FIG. 15. Differential cross section distributions as a function of the invariant mass of the tt pair (left panel) and of total 
transverse energy (right panel) forpp — > tijj +X production at the LHC with ^/s — 7 TeV. The dash-dotted curve corresponds 
to the LO, whereas the solid one to the NLO result. The uncertainty bands depict scale variation. The lower panels display the 
differential K, factor. 



corrections are reduced down to —14%. Again, the LO 
result scales down, this time by 14%, while in the NLO 
case we have an increase of the order of 2%. At LO one 
parton is equivalent to one jet, and the jet algorithm with 
size parameter R is equivalent to the requirement that 
any two partons should be separated by at least ARjj 
distance. A larger value of ARjj means smaller cross 
sections. At NLO, there can be two partons in a jet, and 



jets may have some structure. A larger R means that 
the probability of radiation of a parton out of the area 
with interjet distance R is smaller. This can be trans- 
lated into a reduced number of jets with px below the 
cut threshold and higher total NLO cross section. Con- 
trary to the TeVatron case, this has not been overcome 
here by the simultaneous increase of the ARjj separation 
cut resulting in higher total NLO cross sections. At the 
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FIG. 16. Differential cross section distributions as a function of transverse momentum (left panel) and rapidity (right panel) of 
the 1st hardest jet for pp — > tijj + X production at the LHC with yfs = 7 TeV. The dash-dotted curve corresponds to the LO, 
whereas the solid one to the NLO result. The uncertainty bands depict scale variation. The lower panels display the differential 
K, factor. 
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FIG. 17. Differential cross section distributions as a function of transverse momentum (left panel) and rapidity (right panel) of 
the 2nd hardest jet for pp — > ttjj + X production at the LHC with ^/s = 7 TeV. The dash-dotted curve corresponds to the LO, 
whereas the solid one to the NLO result. The uncertainty bands depict scale variation. The lower panels display the differential 
K, factor. 



LHC, the two hard jets are more likely to be produced 
back-to-back. Therefore, events are mostly concentrated 
around A.Rjj = Tt, which is mildly affected by a change 
in ARjj cut from 0.5 to 1. 

Similarly to the TeVatron case, NLO cross sections for 
two different values of the a ma x parameter are given in 
Table |Xj The independence of the final result on a max is 



obtained again at the permil level. 

In Table IXj predictions for LO and NLO cross sections 
are presented with the default LHC selection, however, 
the transverse momentum cut of the two hard jets is var- 
ied between 50 GeV and 125 GeV. Within a 50 - 100 
GeV range corrections are quite stable, changing the /C 
factor by less than 9%. For 125 GeV cut, we observed 
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FIG. 18. Differential cross section distributions as a function of the invariant mass, rrijj, of the 1st and the 2nd hardest jet and 
Aiijj separation for pp — > tijj + X production at the LHC with y/s = 7 TeV. The dash-dotted curve corresponds to the LO, 
whereas the solid one to the NLO result. The uncertainty bands depict scale variation. The lower panels display the differential 
K, factor. 



a somehow higher rise of NLO QCD corrections up to 
—40%. Compared to the total NLO tt cross section, 

o$%(pp^m = 158.1^; 2 5 jSj Pb, (19) 

we discover that for the smallest pTj cut of 50 GeV the 
tijj events represent only 6% of the total cross section. 
The fraction is decreased to 3%, 1% and 0.6% for 75 GeV, 
100 GeV and 125 GeV respectively. 



2. Differential cross sections 

In the following, we turn our attention to the size of the 
NLO QCD corrections to the differential distributions at 
the LHC with y/1 = 7 TeV. For the LHC, a pp collider, 
the shapes of some distributions are quite different as 
compared to the TeVatron case. Top quark production 
at the LHC, for example, is forward-backward symmet- 
ric in the laboratory frame as a consequence of the sym- 
metric initial state. Thus, Ap B = _4 FB = 0. This can 
be observed in Figure [TU where top and anti-top quark 
rapidity distributions, symmetric around yt, are plotted 
at LO and NLO. Distributions get broadened, moving 
from the TeVatron to the LHC case, as can be expected 
from an increased scattering energy. On the other hand, 
the distribution in the averaged transverse momentum, 
pr t , of the top and anti-top, shown in Figure [121 be- 
comes harder during this transition. Again, the dash- 
dotted curve corresponds to the LO, whereas the solid 
one to the NLO result. At LO the bands correspond to 
the choice of the renormalization and factorization scale 



(i = m t /2 and /i = 2m t . In case of NLO they are de- 
termined out of the following set /i — [mt/2, mt, 2m*]. 
More precisely for each histogram, bands are constructed 
bin-by-bin, by calculating maximal and minimal values 
using the above three scales. Since the scale dependence 
of the NLO cross section at the LHC is a function that 
is not monotonic and the true maximum is somewhere 
close to the central value of the scale, m t , this seems to 
be the only meaningful approach, if one is interested in 
plotting scale uncertainty bands without making a very 
costly scan. Let us emphasize here, that the scale vari- 
ation is, by no means a rigorous way to estimate the 
theoretical uncertainty. At best, it might only give an 
indication of the full uncertainty which is due to the not 
yet calculated higher order corrections. The lower panels 
display the differential K. factor. Angular distributions 
exhibit negative corrections of the order of 30% in the 
central rapidity range and they are decreased to 20% at 
the forward and backward parts of the spectrum. As for 
the averaged transverse momentum distributions of the 
top and anti-top up to 60% distortions are observed. 

In Figure Q31 and Figure [15l differential cross section 
distributions as function of rapidity, y t t, transverse mo- 
mentum, PT tl , and the invariant mass, m t i, of the tt pair 
are given. Furthermore, the total transverse energy of 
the system, Ht, is shown there. The rapidity distribu- 
tion receives negative, but rather stable 25% corrections. 
Unlike pr ti , m-a an d Ht spectrums, where distortions of 
the order of 50%, 35% and 80% respectively, are visible. 

Subsequently, in Figure \W\ and in Figure [T7] differen- 
tial cross section distributions as a function of transverse 
momentum, pi., and rapidity, yj, of the 1st and the 2nd 
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hardest jets are plotted. As anticipated for rapidity dis- 
tributions, in both cases negative and quite stable cor- 
rections of the order of 30% appear. On the other hand, 
for the spectrum, 60% distortions are obtained as 
compared to only 10% in px j2 case. 

Eventually, differential cross section distributions as a 
function of the dijet invariant mass, rrijj, are given in 
Figure [18] Also shown in Figure [18] is the separation 
in ARjj between two jets ordered in hardness. In both 
cases, rather constant —30% corrections can be reached. 
Again, fluctuations visible at the tails of the ARjj dis- 
tribution are caused by limited statistics of the Monte 
Carlo integration. As we can see at the LHC, the two 
hard jets are more likely to be produced back-to-back, 
leading to a more peaked distribution around n. 

Generally, we can say that for a fixed scale fi — fiR — 
= m t at the LHC, the NLO QCD corrections are 
always negative and protean. In case of angular distribu- 
tions, y, AR, and the invariant mass distribution of the 
tt pair, rather steady 20% — 30% corrections are visable. 
In the transverse momentum distributions, on the other 
hand, they introduce up to 60% distortions. In addition, 
in the case of the total transverse energy, Ht, even 80% 
deformations can be observed. 



IV. CONCLUSIONS 

In this paper, we have presented a computation of the 
NLO QCD corrections to tt pair production in associa- 
tion with two hard jets at the TeVatron and the LHC. 
The total cross sections and their scale dependence, as 
well as several differential distributions have been given. 
The NLO QCD corrections reduce the scale uncertainty 
of the total cross sections and of the differential distribu- 
tions compared to LO calculations, which can only pro- 
vide qualitative predictions. 

In case of the TeVatron, the forward-backward asym- 
metry of the top quarks has been calculated for the 
first time. We have found that, with inclusive selec- 
tion cuts, the forward-backward asymmetry amounts to 
4blo = -10-3% at LO and Ap B NLO = -4.6% at 
NLO. 

Furthermore, the impact of the NLO QCD corrections 
on integrated cross sections at the TeVatron is moder- 
ate, of the order of —24%. A the LHC we have obtained 
NLO QCD corrections at the level of -27%. In both 
cases, NLO QCD corrections are reduced substantially, 
if a higher jet resolution parameter R if chosen instead. 
Changing R from R — 0.4 to R — 0.8 in the TeVatron 
case and from R = 0.5 to R — 1.0 in the LHC case de- 
creases NLO QCD corrections down to —16% and —14% 
respectively. A study of the scale dependence of our 
NLO predictions indicates that the residual theoretical 
uncertainty due to higher order corrections is 21% for 
the TeVatron and 15% for the LHC. 

As a further matter, at the TeVatron employing a fixed 
scale fi — /j,r ~ pp = m t , the NLO corrections to differ- 



ential distributions are negative and rather substantial. 
They do not simply rescale the LO shapes, but induce 
distortions up to 60%. Also in the LHC case the NLO 
QCD corrections to the distributions are always negative 
and protean. In the worst case even 80% deformations 
can be observed. 

Finally, our paper demonstrates the utility of the off- 
shell methods and the OPP reduction procedure, as im- 
plemented in the Helac-NLO system, in computing 
NLO QCD corrections for processes of phenomenolog- 
ical interest at hadron colliders. The system consists 
of, CutTools and Helac-ILoop which handle the vir- 
tual corrections and Helac-Dipoles for the real emis- 
sion contributions. For the phase space integration the 
Kaleu package is used and results are cross checked with 
the help of the Phegas phase space generator. More- 
over, OneLOop is deployed for the evaluation of the 
one-loop scalar functions. With this framework we can 
calculate NLO QCD corrections to the remaining multi- 
particle background and signal processes which are rele- 
vant in the ongoing analyzes at the TeVatron and at the 
LHC. Besides, we are prepared to compare NLO QCD 
predictions with already available data. 



V. APPENDIX 

In order to facilitate a comparison to our calculation, 
we provide values of the one-loop virtual corrections to 
the squared matrix elements for one phase space point 
which is given in Table IXIIl Virtual corrections come 
from the interference of the tree level and one-loop am- 
plitudes, summed over all colors and spins, and for rif = 5 
massless quark flavors. As mentioned before we describe 
the running of the strong coupling constant with two-loop 
accuracy, including five active flavors which corresponds 
to a s (m t ) = 0.109403691593462 where m t = 173.3 GeV. 
In Table IXIIII we present numbers for all eight subpro- 
cesses. The quantities presented in Table lXHll correspond 
to the renormalized virtual one-loop correction in the 't 
Hooft-Veltman [28| version of the dimensional regular- 
ization scheme. The remaining singularities in the di- 
mensional regularization parameter d = 4 — 2e arise from 
the virtual soft and collinear singularities in the one-loop 
amplitudes. Additionally coefficients for color and spin 
summed results for the X-operator of the dipole subtrac- 
tion function as defined in [38| are presented. We can 
observe numerical cancellations of infrared divergences 
between the X-operator and the corresponding virtual 
contributions at least up to 10 digits. 
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TABLE XII. The set of momenta for 12 — > tt 3 4 process with the notation (p x ,p y ,p z ,E) and all the components given in Ge V. 
The top quark mass is set to mt = 173.3 GeV. 



Parton 



Py 



Pz 



E 



Pi 
Pi 
Pt 
Pi 

P3 
Pi 






-715.3340594013138 
-24.14783219601748 
21.87826794850173 
717.6036236488295 






-475.1625187999429 
-6.328336607570634 
1000.411563795763 
-518.9207083882493 



2424.746502697520 
-2424.746502697520 
101.1925816377932 
295.5085181344487 
1341.334408905234 
-1738.035508677476 



2424.746502697520 
2424.746502697520 
881.9042263139404 
343.4841188963524 
1673.463460042633 
1950.641200142114 



TABLE XIII. Numerical values of the virtual corrections to the squared matrix elements at the phase space point given in Table 
\XI]\ and = /if = m t for all the partonic subprocesses for the tijj production at a hadron collider. We give the finite parts 
along with the coefficients of the poles in e. Additionally coefficients for color and spin summed results for the T-operator are 
presented. 



SUBPROCESS 



Ve 2 



1A 



l/e° 



99 -> ttgg 
1(e) 



-4.759403785663033x10" 
4.759403785663084x10" 



1.868033206671699x10" 
-1.868033206665652x10" 



-3.727507250164706x10" 



gg —> ttuu 
1(e) 



-1.243597961079720x10" 
1.243597961079740x10" 



4.205081694102065x10" 
-4.205081694097688x10" 



-8.435931194456964x10" 



uu — > ttgg 
1(e) 



-3.408532130149515x10" 
3.408532130149502x10" 



1.296674677886268x10" 
-1.296674677886843x10" 



-2.908105638419785x10" 



uu — > ttuu 
1(e) 



-5.610006218609934x10" 
5.610006218609923x10" 



3.140716653723494x10" 
-3.140716653723457x10" 



-6.264131939313179x10" 



ug — > ttug 
1(e) ' 



-2.178331428688374x10" 
2.178331428688390x10" 



9.512663847206611x10" 
-9.512663847144086 x 10" 



-1.925455596882764X 10 -8 



gu — > ttug 
1(e) ' 



-5.212043686490916x10" 
5.212043686491147x10" 



1.749451737240981x10" 
-1.749451737238111x10" 



-3.143378807927551x10" 



ud ttud 
1(e) 



-5.488979123409352 xlO~ iu 
5.488979123409385 xlO" 10 



3.020380070883476x10" 
-3.020380070883462x10" 



-5.979795310554041 x 10~ 9 



uu — > ttdd 
1(e) 



-2.663317234536143x10" 
2.663317234536076x10" 



1.278242588447108x10" 
-1.278242588446499x10" 



-4.283899547423337X 10" 
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